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04 , Astrophysical plasmas in accretion discs are usually treated in the framework of fluid or MHD 

approaches but there are some situations where these treatments become inadequate and one needs 
to revert to the more fundamental underlying kinetic theory. This occurs when the plasma becomes 
effectively collisionless or weakly-collisional such as, for example, in radiatively inefficient accretion 

■ flows onto black holes. In this paper, we lay down the basics of kinetic theory in these contexts. In 
particular, we formulate the kinetic theory for quasi-stationary collisionless accretion disc plasmas in 
the framework of a Vlasov-Maxwell description, taking the plasma to be non-relativistic, axisymmet- 

■ ric, gravitationally-bound and subject to electromagnetic fields. Quasi-stationary solutions for the 
kinetic distribution functions are constructed which are shown to admit temperature anisotropics. 
The physical implications of the theory are then investigated and the equations of state and angular 

*£h ■ momentum conservation law are discussed. Analysis of the Ampere equation reveals the existence 

of a quasi-stationary kinetic dynamo which gives rise to self-generation of poloidal and azimuthal 
^ ' magnetic fields and operates even in the absence of turbulence and/or instability phenomena. 

I. INTRODUCTION 

Fluids are usually treated as continuous media with their time evolution being determined by a suitable set of fluid 
^ ' equations, even if the underlying fundamental description is in terms of a discrete system of particles, whose dynamics 
is deterministic but characterized by random initial conditions. A continuum description for such systems applies when 
the kinetic description has a continuous phase-space probability density, given by a single-particle kinetic distribution 
function (KDF), which satisfies a suitable kinetic equation. Once the KDF is prescribed, all of the continuum fluid 
moments can be represented in terms of well-defined constitutive equations, determined via appropriate velocity 
moments of the KDF. When binary Coulomb collisions are negligible (collisionless or weakly-collisional plasmas), one 
should return to the kinetic-theory description. In these cases, "stand-alone" fluid or magneto-hydrodynamics (MHD) 
approaches formulated independently of an underlying kinetic theory can usually provide, at best, a partial description 
of the plasma phenomenology. This is because of two possible inconsistencies which may arise. Firstly, the set of fluid 
equations is generally not closed, and so requires the independent prescription of equations of state which may not give 
rise to a self-consistent system. Secondly, in these approaches no account is usually given of microscopic phase-space 
particle dynamics (including single-particle conservation laws) or of phase-space plasma collective phenomena (kinetic 
effects). These issues are naturally addressed within a kinetic treatment, where the fluid fields are determined from 
the underlying self-consistent KDF. This means that both the equations of state and the constitutive equations for 
the fluid fields (see Sections 7-9) then follow uniquely from the microscopic dynamics. 

This paper is concerned with the description of stationary or slowly-time-varying (quasi-stationary) phenomena 
occurring in astrophysical plasmas, focusing particularly on accretion discs around compact objects. For matter 
in the inner parts of such discs, the plasma may be either collisional (as in the majority of cases) or collisionless, 
depending on the circumstances. A notable example of the latter case is provided by radiatively inefficient flows 
(RIAFs) arising in low-density geometrically-thick discs around black holes ([l|). Matter in these systems is thought 
to consist of a two-temperature plasma, with the ion temperature being much higher than the electron one, and the 
Coulomb collision timescale being much longer than the inflow time. Other interesting applications for collisionless 
or weakly-collisional plasmas concern accretion discs around neutron stars and white dwarfs. In the inner regions of 
such discs, the magnetic field of the central object may become dominant and ions and electrons can be collisionally 
decoupled so as to sustain different temperatures. This happens, in particular, if the radiative cooling time-scale of 
the electrons is much shorter than the time-scale for electron-ion collisions so that the electrons and ions are thermally 
decoupled and can have different temperatures (@, Q). 
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Treatments of accretion discs have often been made in terms of a purely fluid-dynamical approach to which was 
added an "anomalous" form of viscosity (i.e. one not due to binary particle collisions), following simple intuitive 
turbulence models such as that giving rise to a-discs. However, it is now almost universally believed that the origin 
of the effective viscosity lies in magnetic phenomena (such as the magneto-rotational instability, MRI) and that the 
medium needs to be treated as a magnetized plasma, when making detailed investigations, rather than as a simple 
unmagnetized neutral fluid. Almost always, these calculations are then performed within the stand-alone continuum 
MHD treatments. However, as indicated above, the underlying plasma kinetic theory needs to be directly considered 
when treating collisionless or weakly-collisional plasmas. An approach of this type has been presented in two recent 
papers (0,01, hereafter referred to as Paper I and Paper II respectively). 

The aim of the present paper is to address in detail the astrophysical aspects of the theory previously developed. 
We focus on accretion discs composed of collisionless plasma which can be subject to strong electro-magnetic (EM) 
fields and which are gravitationally bound, in the sense of being confined by the combined gravitational potential of 
the central compact object and the disc. We deal with plasma which is sufficiently far from the central object so 
that it can be treated non-relativistically, with both kBT/m s c 2 and (i?0/c) 2 being < 0.1, where ks is the Boltzmann 
constant, T is the temperature, m s is the rest-mass of the species of plasma particle labelled with the index s (s = e,i 
for an electron- ion plasma), R is the cylindrical radial distance from the central object (z will be the vertical cylindrical 
coordinate), and is the local angular velocity in the disc. We do not include any interactions with radiation since 
we are not considering any background radiation field and radiation-reaction effects (associated with the self-emission 
of radiation) are negligible for our non-relativistic plasma. We envisage magnetic fields, on the other hand, which can 
in principle be either external (B ea:t ), or self-generated within the disc itself (B se ^). Estimated values coming from 
observations of accretion discs around compact objects include magnetic-field strengths in the range B ~ 10 1 — 10 8 G 
(@i 0)- For these systems, estimates for species temperatures usually lie in the ranges Tj ~ 10 4 — 10 12 A" and 
T e ~ 10 4 — 10 8 A for ions and electrons respectively. Particle densities in different types of accretion disc around 
compact objects span a very wide range of values, but here we focus on the bottom end of it. In particular, for 
discussing the collisionless regime of non-relativistic AD plasmas, the number density is taken to be in the interval 
n s ~ 10 6 — 10 15 cm~ 3 . In the case of hydrogen-ion systems, the Spitzer collision time ra, the Larmor rotation time 
and the Langmuir time r p i are then respectively in the ranges tqi — 10 -2 — 10 19 s, th ~ 10~ 4 — 10 s and 
T p i ~ 1 — 10~ 4 s, while the corresponding estimates for the Debye length Xd, the mean free path X m fp,i and the 
Larmor radius Tl% are respectively Ad — 10~ 2 — 10 6 cm, X m fp,i — 10 6 — 10 31 cm (an extremely large number!) and 
rLi ~ 10~ 4 — 10 7 cro (the lower values corresponding to lower temperature and higher density in the first and second 
cases and to lower temperature and higher magnetic field in the third one) . For a plasma consisting of ion and electron 
species (s = i, e), one can determine the range of values of these characteristic parameters for each species. For any 
phenomena occurring on timescales At and lengthscales AL, satisfying 



these plasmas can be considered as: 

• (#1) Collisionless: due to the inequalities between At and tc s and between AL and A m / P;S . Any plasma 
is effectively collisionless for processes whose timescales and lengthscales are short enough and, under those 
circumstances, one needs to use kinetic theory. For a plasma which is sufficiently diffuse, this can include almost 
all relevant processes. 

• (#2) Characterized by a mean-field EM interaction: due to the inequalities between At and r ps and between 
AL and Xd, charged particles of the accretion disc interact with the others only via a continuum mean EM 



• (#3) Quasi-neutral: due to the inequality between AL and Ad, the plasma can be taken as being quasi- neutral 
on the lengthscale AL. 

When conditions #1 and #2 are satisfied, the medium is referred to as a Vlasov- Maxwell plasma and kinetic theory 
needs to be used. It is this which we will be considering in the following. The plasma is treated as an ensemble 
of particle species, each being described by a KDF f s (y,t) where y is the particle state vector and t is the time, 
with y = (r, v) where r is the position and v is the velocity. The species KDFs satisfy the Vlasov kinetic equation 
-iffs(y,t) = 0, with the velocity moments of the KDFs determining the source of the EM self-field {E se ^,B se ^}, 



Tps,T Ls < At < T Cs , 

A fl « AL« X mfPr 



(1) 
(2) 



field. 
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identified with the plasma charge and current density {p(r, t), J(r, t)}: 



P {v,t) = Y, z * e rf Vs(y^)> 

s ^ 

J(r,t) = ]Tz e e / d 3 vvf s (y,t) 



(3) 



(4) 



We work here within the framework of perturbative kinetic theory, extending the approach developed in Papers I 
and II. This approach allows construction of analytical solutions for the KDF for slowly-time-varying axisymmetric 
gravitationally-bound systems (referred to as quasi- stationary KDFs). The treatment is developed for both strongly 
and weakly-magnetized plasmas, distinguishing between magnetic field configurations with closed and open magnetic 
surfaces ([8l4ll|V In the present paper, we are focusing on equilibrium or quasi-stationary configurations as a pre- 
liminary to subsequently studying perturbations around these solutions. Note that what is meant here by the term 
"equilibrium" is in general a quasi-stationary KDF, expressed in terms of the relevant first integrals of motion and adi- 
abatic invariants (see Papers I and II and |10l - [l3| ). which can also include a non- vanishing stationary radial accretion 
flow ( 14]). One of our main results is the demonstration that for strongly- magnetized collisionless plasmas, kinetic 
theory gives the possibility of having quasi-stationary accretion flows even in the absence of any additional forms 
of effective viscosity, such as those arising from turbulence phenomena (which lie beyond the equilibrium solutions). 
Within the framework of kinetic treatment, quasi-stationary accretion flows can occur in collisionless AD plasmas 
as part of the equilibrium configuration, with the role of viscous stresses being played by the anisotropic pressure 
tensor generated by phase-space anisotropies. We note that a somewhat similar magnetically-driven non-turbulent 
mechanism for driving accretion has been pointed out before ([HI), based on a purely fluid treatment. While that has 
some similarity with what is being described here, in the sense that magnetic phenomena in stationary configurations 
are giving rise to the redistribution of angular momentum required for having the accretion flow, it is very different 
in other respects. 

We are concerned here particularly with investigating the role of specifically kinetic effects which do not appear in a 
fluid treatment. These can arise due to individual-particle dynamics (e.g. finite Larmor-radius effects and ones arising 
from microscopic conservation laws), as well as due to statistical properties of the equilibrium KDF (e.g. temperature 
anisotropy and non- uniform fluid fields). Their specific influence is found to depend critically on the topology of the 
magnetic surfaces, the statistical properties of the plasma and the strength of the EM field. Instabilities and turbulent 
phenomena (see for example [161 Il7||) are not considered here. 

A number of issues arise for accretion-disc plasmas of this type, including ones related to: 

1) Existence of kinetic equilibria for both strongly and weakly-magnetized plasmas; 

2) Dynamo effects occurring in quiescent accretion-disc plasmas, which can explain the self-generation of both 
azimuthal and poloidal magnetic fields; 

3) Collisionless non-turbulent quasi-stationary accretion processes; 

4) Closure conditions yielding a finite set of fluid equations for the relevant fluid fields; 

5) Equations of state for the pressure tensor components. 

Investigation of these issues is relevant for correctly understanding the equilibrium and dynamical properties of 
collisionless accretion-disc plasmas. 

The paper is organized as follows. In Section 2 we provide a classification of accretion-disc plasmas together with the 
basic assumptions and definitions being used. In Section 3, the first integrals and adiabatic invariants of the system are 
derived, and their physical meaning is discussed. Section 4 deals with the construction of the quasi-stationary KDF 
for strongly-magnetized plasmas with open magnetic surfaces. Similar calculations are then presented in Section 5 for 
strongly-magnetized plasmas with closed nested magnetic surfaces, and in Section 6 for weakly-magnetized plasmas. 
In Section 7 the number density and flow velocity of each species are computed for all of the configurations considered 
and their physical properties are analyzed, stressing their connection with accretion-disc dynamics. In Section 8 the 
pressure tensor for each species is explicitly calculated. A separate discussion is provided concerning the determination 
of equations of state for the pressure tensor components. Section 9 deals with the calculation of the fluid angular 
momentum conservation law, while Section 10 deals with the kinetic accretion process which does not depend on the 
presence of perturbative processes or an explicit viscosity. Section 11 contains a demonstration of the existence of a 
kinetic dynamo mechanism and presents equations giving the poloidal and toroidal components of the magnetic field. 
Finally, Section 12 contains a summary of the main results and closing remarks. In the technical parts of this paper, 
we will often abbreviate "accretion-disc plasma" to "AD plasma" for conciseness. 
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II. BASIC PHYSICAL ASSUMPTIONS 



In the following we will distinguish between strongly-magnetized, intermediate and weakly-magnetized accretion- 
disc plasmas, the distinction being made on the basis of asymptotic conditions expressed in terms of suitable species- 
dependent dimensionless physical parameters. These parameters will also be used for characterizing the EM fields 
and for constructing perturbative kinetic solutions for the quasi-stationary KDFs describing AD plasmas. Three 
parameters are required: 

1) The first parameter is Sm.s = where s = i,e again denotes the species index. Here r^ s = v± t hs/^cs is the 
species average Larmor radius, with v_i_ths = {ksT^s/Ms} 1 ^ 2 denoting the species thermal velocity perpendicular to 
the magnetic field and O cs = Z s eB /M s c being the species Larmor frequency. Note, that we are here always measuring 
temperatures in degrees Kelvin, with ks denoting the Boltzmann constant. L is the characteristic length-scale of 
the spatial inhomogeneities of the EM field, defined as L ~ Lb ~ Le, where Lb and Le arc the characteristic 
magnitudes of the gradients of the absolute values of the magnetic field B (x, t) and the electric field E (x, t) , defined 



as -Y- 



1, 3 j and j-^ = max j | In E , i = 1, 3j, where the vector x denotes x = (R, z). From 



Sm.s it is possible to define a unique parameter sm = rnax{£M, s , s = i, e}. 



2) The second parameter is defined as e s 



M s Rv^ 



i.e. it is the ratio between the toroidal angular 



momentum of the particle L vs = M s Rv v and the magnetic contribution to its toroidal canonical momentum p vs , 
with p vs = L vs + ^^ip- Here v v = v • e v , with e v being the unit vector along the azimuthal direction tp, while ip 
denotes the flux function of the poloidal magnetic field (see definition below). 

3) The third parameter is <r s = 



M s 2 




k B T 3 






Z s eQ e s " 



which represents the ratio between the kinetic energy of the 

particle ^-v 2 and its potential energy. The latter is represented in terms of the effective potential <&^, which, in 
turn, is related to the electrostatic and gravitational potentials (see the next section for a rigorous definition), while 
for thermal plasmas one has Mf-v 2 <~ ksT s . Related to a s , we also define a = max (cr s ,s — i,e). The AD plasma is 
said to be gravitationally bound if a <C 1 , and we are only considering such cases in the following. 
We classify these AD plasmas to be: 

1) Strongly-magnetized if the two conditions Sm.s <C 1 and e s 1 are satisfied. 

2) Intermediately-magnetized if £m,s <C 1 but e s ~ 1. 

3) Weakly-magnetized if em,s ~ 1 and e s > 1. 

In both the strongly and intermediately magnetized cases, charged particles spiral tightly around magnetic field 
lines (so that gyrokinetic theory applies, see below), whereas they do not do this in weakly- magnetized plasmas. In 
this sense, the behaviour of strongly and intermediately magnetized plasmas is similar but, on the other hand, they 
are importantly different in that the magnetic flux ip completely dominates the particle canonical momentum for 
strongly-magnetized plasmas, so that the contribution due to the particle angular momentum L vs is negligible in the 
corresponding conservation law, whereas the contributions are comparable for intermediately-magnetized plasmas. 

To give an indication of the circumstances in which conditions 1) and 2) above apply, we consider two situations 
representative of stellar mass and galactic-centre mass black holes. Taking mass M* = 10M© (giving a Schwarzschild 
radius Rs c h ~ 30km) as representative of a stellar-mass black hole, we focus on AD plasma located at R ~ 10 Rsch 
from the central object, with ion temperature Tj ~ lO 8 ^, and characteristic length-scale L ~ 0.1R = Rs c h for the 
EM field. Then taking sm,i ^ 10~ 3 as representative of Sm,i < 1, we get B > 10 G which is at the lower end 
of the range of magnetic field strengths given above, and so the condition ema -C 1 can easily be realized. For a 
galactic-centre black hole with mass = 1O 8 M , the equivalent estimate, for the same radial distance in terms of 
Schwarzschild radii, gives B > 10~ 6 G. Considering now the condition e s <C 1, we make similar choices as above, plus 
taking the representative azimuthal particle velocity v v as being between the ion thermal velocity and the Keplerian 
velocity, and the magnetic field strength as being related to the magnetic flux ip by B ~ j^. Then requiring, for 
example, < 10~ 2 gives the (very approximate) bounds of B > 10 8 G for the stellar-mass case and B > 10 G for the 
galactic-centre one. The bound for the stellar-mass case is at the upper end of the range for magnetic field strengths 
given above, while that for the galactic-centre one is far more modest (and we note that this is the more interesting 
context in practice for the appearance of RIAF-type accretion flows). However, concerning these estimates, one should 
bear in mind that, in comparison with the case for Sm.s, the bound on B coming from e s is much more sensitive to 
the species being considered (electrons or ions) as well as to the value of the particle angular momentum and to the 
approximate order-of-magnitude relationship used for estimating B from ip. Depending on the circumstances being 
considered, the condition e s <C 1 may or may not be satisfied, and so one should investigate both possibilities. 

From the above considerations, it follows that both strongly and intermediately-magnetized plasmas may be found 
in the inner regions of accretion discs around compact objects. Weakly-magnetized plasmas would be located further 
out, in the outer regions of discs, where lower temperatures and weaker magnetic fields are expected. We will focus here 
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FIG. 1: Schematic view of the topology of the magnetic surfaces. 



on studying the configurations for strongly and weakly-magnetized collisionless AD plasmas, leaving the corresponding 
investigation of intermediately-magnetized plasmas to a separate study. 

Ignoring possible weakly-dissipative effects (Coulomb collisions and turbulence), we will assume that the KDF and 
the EM fields associated with the plasma obey the system of Vlasov-Maxwell equations, with Maxwell's equations 
being considered in the quasi-static approximation. For definiteness, we will consider here a plasma consisting of two 
species of charged particles: one species of ions (i) and one of electrons (e). 

Following the treatment presented in Papers I and II, we are taking the AD plasma to be: a) non-relativistic, 
in the sense that it has non-relativistic species flow velocities and fcsT/m s c 2 , that the gravitational field can be 
treated within the classical Newtonian theory, and that the non-relativistic Vlasov kinetic equation can be used as the 
dynamical equation for the KDF; b) collisionless, so that the mean free path of the plasma particles is much longer 
than the largest relevant characteristic scale length of the plasma; c) axisymmetric, so that the relevant dynamical 
variables characterizing the plasma (e.g., the fluid fields) are independent of the azimuthal angle tp, when referred to 
a set of cylindrical coordinates (R, ip, z); d) acted on by both gravitational and EM fields. 

We focus here on solutions for the equilibrium magnetic field B which admit a family of locally-nested axisymmetric 
toroidal magnetic surfaces ( [1, H[ ) . We recall that a magnetic surface is defined as a surface on which the poloidal 
magnetic flux ip is constant, and that the condition V?/> • B = is then identically satisfied on each magnetic surface. 
A schematic view of nested magnetic surfaces is shown in Figure 1. For the equilibrium configurations which we are 
considering here (prior to any subsequent perturbation) it is reasonable to think that the magnetic field would be 
rather ordered, at least on a local scale, so that having nested magnetic surfaces is likely on that local scale. We 
distinguish between the cases of locally closed magnetic surfaces, discussed in Paper I, and locally open magnetic 
surfaces, discussed in Paper II. Note that here the meaning of open and closed surfaces has to be interpreted with 
reference to the local domain occupied by the AD plasma. See Figure 1 for a schematic comparison of the two 
topologies. For both of the configurations, a set of magnetic coordinates (if!,(p,'&) can be defined locally, where d is 
a curvilinear angle-like coordinate on the magnetic surfaces f/'( x ) = const. Each relevant physical quantity G(x, t) 
can then be conveniently expressed either in terms of the cylindrical coordinates or as a function of the magnetic 
coordinates, i.e. G(x,t) = G (V>, i) , where the ip dependence has been suppressed due to the axisymmetry. 

For the configuration of closed nested magnetic surfaces, we will assume "small inverse aspect ratio ordering" to 
hold, as expressed by the requirement that < 5 <C 1. Here the dimensionless quantity S = ^-jf^ is referred to 
as the inverse aspect ratio parameter, where Rq is the radial distance from the vertical axis to the centre of the 
nested magnetic surfaces and r max is the average cross-sectional poloidal radius of the largest closed toroidal magnetic 
surface. This ordering is consistent with the results presented by [8l-fll| and with the assumption of nested and closed 
magnetic surfaces which are assumed to be localized in space (see also further discussion in Paper I). 

In the following we will consider AD plasmas which are characterized by slowly time-varying phenomena. In 
particular, the EM field is taken to be given by an analytic function of the form 



with k being an integer > 1 and A = min (em-, &)■ This time dependence is connected with either external sources or 



[E(x,A fc t),B(x,A fc i)] , 



(5) 
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boundary conditions for the KDF. In particular, we will assume that the magnetic field is of the form 

B = V x A = B sei/ (x, X k t) + B ext (x, X k t), (6) 

where B seZ ^ and B eK * denote the self-generated magnetic field produced by the AD plasma and a finite external 
magnetic field produced by the central object (in the case of neutron stars or white dwarfs). For greater generality, 
we will not prescribe any relative orderings between the various components of the total magnetic field, which are 
taken to be of the form 

B self = J(x, X k t)\7p + W> p (x, \ k t) x Vtp, (7) 
B ei ' = Vfc(x,A' : f)xV^. (8) 

In particular, here Bj = 7(x, X k t)\7ip and Bp e Vip p (x, X k t) x \7ip are the toroidal and poloidal components of 
the self-field, while the external magnetic field R ext is defined in terms of the vacuum potential ^>£>(x, X k t). As a 
consequence, the magnetic field can also be written in the equivalent form 

B = I(x, \ k t)\7tp + W>(x, X k t) x V</5, (9) 

where the function tp(x,X k t) is defined as ip(x,X k t) = ip p (x,X k t) + V'd(x, X k t), with k > 1 and (ip,ip,$) defining a 
set of local magnetic coordinates (as implied by the equation B • Vip — which is identically satisfied). Finally, 
it is also assumed that the charged particles of the plasma are subject to the action of effective EM potentials 
{$^/(x, X k t), A(x, X k t)} , where $ e / f {x, X k t) is given by 

M 

*«"(x, X k t) = $(x, X k t) + -^<Mx, A fc t), (10) 

Zj s e 

with $!^(x, X k t), <&(x, X k t) and $G( x )A fe i) denoting the effective electrostatic potential and the electrostatic and 

gravitational potentials. In particular, the gravitational potential $g is expressed as $g = <&^ xt ^ + $^ ei ^- ) , where 

<f>Q Xt ^ and denote the gravitational potentials due to the external sources (i.e., primarily the central object) 

and the same accretion disc respectively. The latter can be computed from the Poisson equation 

V 2 $^ (/) (x, X k t) = 4nG NPm (x, X k t) , (11) 

where Gjv is the Newton gravitational constant and the source term p m (x, X k t) describes the distribution of matter 
of the AD plasma. In the following it is assumed that the gravitational potential $g is a known quantity. 



III. FIRST INTEGRALS AND ADIABATIC INVARIANTS FOR ACCRETION-DISC PLASMAS 

In this section we first define relevant kinetic conserved quantities, namely dynamical variables depending on the 
states of the individual charged particles of the plasma. We refer here in particular to collisionless AD plasmas 
which satisfy the assumptions introduced in the previous section. Then, on the basis of this, we will investigate the 
basic physical implications of the conservation laws for the qualitative properties of single-particle dynamics. The 
determination of dynamical invariants (i.e., first integrals of motion and adiabatic invariants) is a basic requirement 
for the construction of equilibrium solutions for the KDF (kinetic equilibria). We recall here that an adiabatic 
invariant P of order n with respect to A is a quantity which is conserved only asymptotically, i.e., in the sense that 
jp-^lnP = + 0(A n+1 ), where n > is a suitable integer. In other words, adiabatic invariants are dynamical 
variables which change slowly on the time-scale of the Larmor rotation time. These invariants can be derived from 
the Lagrangian formulation of the single particle dynamics and the corresponding gyrokinetic theory (see subsection 
below), which can be obtained by means of a suitable asymptotic expansion in terms of the dimensionless parameter 
Em,s as long as this is <C 1. Because of this requirement, and in view of the previous classification, gyrokinetic theory 
can only be formulated for strongly-magnetized and intermediately-magnetized AD plasmas. 

We first consider the treatment of the first integrals and adiabatic invariants which are conserved for both strongly- 
magnetized and weakly-magnetized AD plasmas. By assumption, the only first integral of motion is the toroidal 
canonical momentum p vs conjugate to the ignorable azimuthal angle ip: 

p vs = M s Rv • e v + — tp = —^-ip* s . (12) 
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The total particle energy 



M 

E s = ^v 2 +Z a e$ e /f(x,\ n t), (13) 



is instead considered, by construction, to be an adiabatic invariant of order n, with n > 1. 

A. Guiding-centre adiabatic invariants for strongly-magnetized plasmas 

Additional adiabatic invariants can be determined for strongly-magnetized plasmas based on gyrokinctic (GK) 
theory (see Paper II, [18-20]). We recall that gyrokinctic theory provides a convenient formulation for single charged- 
particle dynamics in the presence of EM fields which are strong enough so that the motion of the particle can be 
well-represented as spiraling around a single field line, following an imaginary guiding centre which moves along 
the field line. Gyrokinetic theory clearly displays the characteristic symmetry of the spiral motion, allowing further 
conservation laws to be derived in addition to those mentioned previously. We refer to Paper II for the mathematical 
details concerning the formulation of non-relativistic gyrokinetic theory in the presence of both EM and gravitational 
fields. Here, for the sake of completeness, we recall the basic notation for this and the main conclusions which can be 
inferred from it. In the following, we will use a prime " ' " to denote a dynamical variable defined at the guiding- centre 
position r' (or x' = (R' , z') in axisymmetry) . The single particle velocity is decomposed as 

v = 7i'b' + w , + V^ // , (14) 

where u'h' and w' denote respectively the parallel and perpendicular components of the guiding-centre velocity with 
respect to the magnetic field direction, while the effective velocity V e ff is defined as 

V; // (x > ek*) = -J 7 E'/"xb', (15) 

with b' = b(x',e| f t), b(x,e| f f)= B(x,ey)/S(x, sy) being the unit vector of the local magnetic field and where 
we have identified A = Em- It then follows immediately that an adiabatic invariant is provided by the guiding- centre 
canonical momentum p' . Correct to 0{e\^), with fc > 1, this is given by 



, _ m s ( cvy.y^^ Z s e 

P** = W\ Ul B> + — (16) 



A further adiabatic invariant is the magnetic moment m' s , which is proportional to the canonical momentum p<j,i s 
conjugate to the gyrophase angle 4>' , which can be determined in principle to arbitrary order in em f[2l|). In particular, 

the leading-order approximation yields m' s = fx' s = . Finally, by construction, the guiding- centre Hamiltonian 

T-L' s is also an adiabatic invariant. Accurate to order n this is given by 

%'P ee m' s B' + (u>b> + V> ff ) 2 + Z s e^f f - (17) 

When the above assumptions hold, the invariants determined in this way necessarily exist for arbitrary initial condi- 
tions. 



B. Physical implications of the conservation laws 



We now discuss the physical meaning of the conservation laws introduced here and their implications for particle 
dynamics in magnetized accretion discs. 

Consider first the conservation of the toroidal canonical momentum (fl"2j) . For a charged particle it follows that this is 
the sum of two terms: the particle angular momentum L vs = M s Rv v and a magnetic contribution ^rip- It follows that 
the angular momentum by itself is generally not conserved. As a consequence, the canonical momentum conservation 
law allows for the existence of radial particle motion inside a disc. In fact, since in AD plasmas the magnetic flux 
function ip is necessarily spatially non- homogeneous, a moving particle must change its angular momentum L vs 
while fulfilling the constraint ip<, s — const., namely staying on a V'*s — surface. Depending on the geometry of the 
magnetic surfaces, such particle motion may correspond to either a vertical or radial velocity towards regions of 
higher or lower magnetic flux. Since the magnetic contribution to ip# s depends on the sign of the charges, single ions 
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and electrons exhibit motions in different directions while keeping tj)^ s constant. This feature is different from the 
situation for neutral particles, for which the angular momentum itself is conserved. Because of the presence of plasma 
boundaries, this can lead to the self-generation of quasi-stationary electric fields in the accretion disc as a result of 
charge separation. 

We next focus on the conservation of the guiding-centre Hamiltonian (|17p and the magnetic moment f/ s (to leading- 
order approximation). These can be combined to represent the parallel velocity u' as 



Therefore v! is a local function of the guiding-centre position vector x' and, due to axisymmetry, of the corresponding 
flux coordinates {if}' ',*&')■ The above relationship is the basis of particle trapping phenomena, corresponding to the 
existence of allowed and forbidden regions of configuration space for the motion of charged particles. In fact, since 
u' is only defined in the subset of the configuration space spanned by (?//,"$') where the argument of the square 
root is non-negative, it follows from Eq.(fT8|) that particles must undergo spatial reflections when v! — 0. The points 
of the configuration space where this occurs are the so-called mirror points and the occurrence of such points may 
generate various kinetic phenomena. In particular, particles can in principle experience zero, one or two reflections 
corresponding respectively to passing particles (PPs), bouncing particles (BPs) and trapped particles (TPs). In the 
present case, since the right hand side of Eq. (TT8")) depends on the magnitude of the magnetic field (£>'), the effective 
potential energy (Z s e$'gff) and the centrifugal potential we will refer to the TP case as gravitational EM 

trapping. 

Finally, an important qualitative property of collisionless magnetized plasmas follows from the conservation of the 
magnetic moment fi' s . The expression for this relates the magnitude of the perpendicular velocity w' to that of the local 
magnetic field B' . Conservation of the adiabatic invariant /i' s implies that when a charge is subject to a non- uniform or 
a non-stationary magnetic field, its kinetic energy of perpendicular motion, AI S ^-, must change accordingly so as to 
keep fj,' s constant. On the other hand, particles moving on "0» s — surfaces generally necessarily experience a non-uniform 
magnetic field B(x, e\jt). It can be shown that this property implies also the phenomenon of having a non-isotropic 
kinetic temperature (i.e. there being different effective temperatures parallel and perpendicular to the local direction 
of the magnetic field). From the statistical point of view of kinetic theory, this temperature anisotropy corresponds 
to an anisotropy in the kinetic energy of random motion of particles subject to the magnetic field. Such a feature is a 
characteristic kinetic phenomenon arising in magnetized collisionless plasmas. This physical mechanism operating at 
the level of single particle dynamics has important consequences also for the macroscopic properties of such plasmas. 
As we will see, conservation of fi' s allows the effects of temperature anisotropy to be included consistently in the 
quasi-stationary solution for the KDF, and for its physical implications for the dynamics of the corresponding fluid 
system to be inferred. Another candidate source of temperature anisotropy is radiation emission (cyclotron radiation) 
due to Larmor rotation in the presence of a strong magnetic field. The signature of this is the simultaneous occurrence 
of radiation emission corresponding to the Larmor frequencies of the different plasma species. 

IV. THE QUASI-STATIONARY KDF FOR STRONGLY-MAGNETIZED PLASMAS: THE CASE OF 

OPEN SURFACES 

In this section and the following one, we point out the most relevant physical aspects of the kinetic treatment of 
quasi-stationary strongly-magnetized collisionless AD plasmas in the presence of magnetic field configurations with 
nested magnetic surfaces. In this section, we consider the case of magnetic surfaces which are open in the domain of 
the plasma. The case of closed magnetic surfaces is treated in the following section. Following the treatment presented 
in Paper II, here we want to emphasize the physical aspects of the theory and the role of the kinetic approach adopted 
here. For greater generality, it is assumed that each species in the plasma is associated with a set of sub-species (the 
PPs, BPs and TPs mentioned above), each one having a different KDF. The existence of temperature anisotropy is 
allowed for all of the species, involving the introduction of different temperatures parallel and perpendicular to the 
local direction of the magnetic field. It is also assumed that a non- vanishing species-dependent poloidal flow velocity 
can exist, related to possible inward or outward matter flows in the disc (see also the discussion below). Under these 
assumptions, as pointed out in Papers I and II, a solution for the quasi-stationary KDF can be obtained, which in 
the following is denoted as /* s . This is expressed in terms of the integrals of motion and the adiabatic invariants 
identified in the previous section. It has the form 




(18) 




(19) 
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where k > 1 and A is to be identified with sm- Here, by construction, f* s is only defined in the subset of the phase- 
space where the adiabatic invariants are defined, while the variables in the round brackets (E s , ip* s ) will be involved 
in the perturbative expansions to be defined below. It follows that /* s is suitable for describing passing, bouncing and 
trapped particles. Hence, by construction, the KDF is itself an adiabatic invariant, and is therefore an asymptotic 
solution of the Vlasov equation. Due to the arbitrariness of the definition of the KDF, for each plasma sub-species it 
is always possible to identify it with a superposition of Gaussian distributions or, more generally, suitably generalized 
Gaussian distributions. However, for a collisionless plasma, each of these functions must actually itself be a quasi- 
stationary solution and so, due to the requirement (fTT?|) . / ts can always be prescribed to be asymptotically "close" to a 

local bi-Maxwellian KDF. Regarding this, in Paper II it was proved that /* s can be identified with a properly-defined 
Generalized bi-Maxwellian KDF with parallel velocity perturbations. 

In this paper we are mainly interested in the astrophysical applications of the kinetic analysis and so we omit here 
all of the mathematical details of the derivation, referring to Paper II for an exhaustive discussion. It is sufficient to 
mention here that a characteristic feature of the quasi-stationary KDF f* s is that it contains implicit dependences 
in terms of the single particle velocities. As shown in Paper II, these dependences can be made explicit for strongly 
magnetized and gravitationally bound plasmas so as to allow an asymptotic analytical treatment of the velocity 
moments. 

An interesting feature of the theory is the double Taylor-expansion which is performed on /* s to reach this goal, 
achieving a systematic solution method for the Vlasov equation. More precisely, this is done in terms of the two 
dimensionless parameters e s and cr s , requiring that the inequalities e s <C 1 and o~ s <C 1 are satisfied when the 

1/2 

particle velocity is taken to be of the order |v| < v± t hs or |v| < uny^, where vj_ t h s = {^B^J-s/Ms} an d V\\ths = 

{fes2]|s/ Ms} 1 ^ 2 denote respectively the parallel and perpendicular thermal velocities. For greater generality, e s and 
o~ s are here treated in the perturbative expansions as being infinitesimals of the same order. No expansion is performed 
in the ratio The Taylor expansion of /* s with respect to the two parameters can be formally carried out at any 

order in e s and a s . In particular, to leading-order this is done by setting -0*s — iji + O (e s ) and E s = Z s e<fr e Jf + O (cr s ). 
It is then straightforward to prove that the following relation holds to leading-order for the quasi-stationary KDF: 

U a = f s [l + h 1 ^ + h 2 Ds ] , (20) 

where the leading-order distribution f s is of the form f s = f s [E s ,ip^ s ,p' ,m' s ,(ip,'d) ,X k t). Here h 1 Ds and h 2 Ds 

represent the so-called FLR-diamagnetic and energy- correction parts of /* s which are polynomial functions of the 
particle velocity, while f s can be written as 

fs 



x exp < 



(27rk B /M s f 2 (7|| S ) 1/2 T ±S 



M, ( v — V, — Uuh 1 



2k B T l[s s k B A Ts 



(21) 



which we refer to as the bi-Maxwellian KDF with parallel velocity perturbations. Again we stress that f s is only defined 
in the subset of phase-space where the parallel velocity \u'\ is a real function. Here we note that the form of f s has 
been prescribed in order to allow the existence of: 

1) A finite azimuthal flow velocity V s = e v RQ s , with tt s being a suitable rotational frequency (a toroidal angular 
velocity). 

2) Finite parallel flows (with respect to the local magnetic field direction) associated with U[, s = -gr^s- These can 
include inward or outward radial flows of matter. The guiding-centre parallel flow velocity U!, s is uniquely prescribed 
in terms of £ s , with £ s being a suitable frequency (see Paper II). 

3) A finite toroidal magnetic field, which is related to UL. In fact UL is non- vanishing only for magnetic configu- 
rations in which the toroidal magnetic field is non-zero. 

4) A finite temperature anisotropy identified with kB \ T = k B T ± fegV ' 

5) A non-uniform effective number density defined as 

n s = n s exp , (22) 
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with rj s denoting the pseudo-density. Here the function X s is prescribed in such a way as to take into account the 
effects of the electrostatic and gravitational energy of the particle, the centripetal potential and azimuthal and parallel 
flows (the precise definition of this is given in Paper II). 

6) Separate treatments of species and sub-species contributions, for which the previous asymptotic orderings are 
assumed to hold. In fact, for the different populations the analytical expansion can lead to different contributions 
for the terms appearing in the diamagnetic and energy-correction parts, depending on the relative magnitudes of the 
parameters e s and o~ s . On the other hand, because of the double expansion and the energy dependence, the asymptotic 
solution for the two species can hold also in different spatial domains. 

Finally we note that the expansion given in Eq. (l20]) shows that, in general, a bi-Maxwellian KDF cannot be an 
exact stationary solution of the Vlasov equation. Instead, the actual (asymptotic) equilibrium is necessarily described 

by the quasi-stationary KDF /* s . By construction this is asymptotically close to a bi-Maxwellian when the expansion 
pOjl holds. 

There are two important physical implications which follow from the quasi-stationary solution for the KDF (|20[) . 
The first one concerns the existence of kinetic constraints, namely functional dependences which need to be imposed 
on the quasi-stationary KDF in order to guarantee that this is an adiabatic invariant of the prescribed order. To 
outline this point, consider the set of functions 

a s = (p a ,a;,T h ,n a ,z a ), (23) 

which we will refer to as structure functions (see also Paper II). In particular, here 

Ps = t~tf; j (24) 
U ' (25) 



which depend on the pseudo-density, the magnitude of the guiding-centre magnetic field and the parallel and per- 
pendicular temperatures. It is important to point out that the kinetic constraints actually prescribe well-defined 
functional dependences for the structure functions, imposing for them the form 

A, = A s Z s e$ e /f) + O (e.) + O {tr.) . (26) 

The effective potential $ e Jf is generally a function of the form <& e Jf = <I>^(x, e^t), with x = (R, z), while neither the 
gravitational potential nor the electrostatic potential are expected to be functions only of ip. Therefore, in magnetic 
coordinates, the structure functions are of the general form A s = A s (ip, £ k M i) . Hence, the functional forms of the 
leading-order effective number density, the parallel and azimuthal flow velocities and the temperatures carried by the 
bi-Maxwellian KDF, are uniquely determined in terms of tp and As a consequence, the azimuthal angular velocity 
is of the general form f2 s = fi s (ip,'d,e\ I t). We stress that in the customary treatment of collisionless AD plasmas 
based on ideal-MHD these constraints are missing. Instead, they follow in a natural way from kinetic theory. By 
adopting a kinetic treatment it is possible to prescribe the correct form for the fluid fields, as required by the presence 
of kinetic constraints. 

A further important consequence of the kinetic constraints is the relationship between the magnitude of the tem- 
perature anisotropy and the guiding-centre magnetic field at two different spatial locations. In fact, the quantity 
-^p- in the KDF is necessarily an adiabatic invariant. To leading-order in the GK expansion, this implies that the 
asymptotic equation 

[AtJ 2 „ [B} 2 



(27) 



must hold identically for any two arbitrary positions "1" and "2", with ([Ajjj , [B]^) and ([AtJ 2 , [B] 2 ) denoting the 
temperature anisotropy and the magnitude of the magnetic field at these positions respectively. 

The second striking aspect of the kinetic treatment concerns the diamagnetic and energy-correction contributions 
h\) s and h 2 Ds of /* s . These carry the contributions from the expansions of the particle toroidal canonical momentum 
and particle total energy respectively. The perturbative-correction to the KDF is a polynomial function of the particle 
velocity which depends linearly on the so-called effective thermodynamic forces. The latter are here denoted as Ai S 
and C'i s , with i — 1,5. In analogy with classical thermodynamics, it is natural to identify them with the gradients 
of the structure functions A s . Hence, in the present case they are associated with partial derivatives taken with 
respect to the magnetic flux tp and the effective potential <& e /f (9)- There are the following definitions: A\ s = 9 ^ s , 



11 



^2S - Q$ , _ g^ , A 4s _ -g^", A 5s = —g^JJ- and O is = ^e//' , ^2S = g^J f , ^3S = g^J } , ^ 4s = ^77^ 

° 5s = aW 7 ' 

The diamagnctic and energy-correction effects carried by h Da and h 2 Ds cannot be ignored: the construction of 
kinetic equilibria cannot be achieved without them, as pointed out in Papers I and II. From a physical point of view, 
the perturbative contribution to the KDF determines first-order corrections to the fluid moments of /* s produced 
by finite Larmor-radius (FLR) effects. These carry the contributions of all of the thermodynamic forces which can 
arise in collisionless plasmas characterized by non-uniform differential rotation, density and temperature gradients 
and temperature anisotropy. 



V. THE QUASI-STATIONARY KDF FOR STRONGLY-MAGNETIZED PLASMAS: THE CASE OF 

CLOSED SURFACES 

In this section we present the kinetic solution for strongly-magnetized plasmas in the case in which the equilibrium 
magnetic field admits locally a family of closed and nested magnetic surfaces. This is the configuration considered in 
Paper I; the geometry is illustrated schematically in Fig.l. 

The quasi-stationary KDF for collisionless plasmas with closed magnetic surfaces can be found as a particular 
asymptotic limit of the general solution /* s holding for open surfaces. Besides considering closed surfaces, this is 
obtained by imposing the requirement of small inverse aspect ratio ordering. As a side assumption, we must now 
also impose vanishing of the velocity perturbation Ul, which requires setting £ s to zero. In fact, the case of closed 
surfaces corresponds to plasma magnetic self-confinement in which no local net radial flow can take place, the latter 
being associated with U'^ s (see Section 10). Note however that a general situation can include both locally-closed and 

open magnetic surfaces, as shown in Fig.l. In this context, the quasi-stationary KDF /* s is reduced to a Generalized 
bi-Maxwellian KDF of the form 

Us = Us {E s ,i>* s ,m' s , (ip, s ),X k t) , (28) 

with k > 1 and A being identified with em, while the perturbative expansion is applied only to the variable (*p* s ). When 
our set of assumptions holds, the Taylor expansion of f* s can be performed only with respect to the dimensionless 
parameter e s , while again no expansion is performed in the ratio —^5. Also, it is possible to prove that in the present 
case the perturbative contributions with respect to the expansion in a s all become negligible. More precisely, correct 
to first-order in e s , the asymptotic expansion of /* s gives: 

Us = f s [1 + h 1 ^] , (29) 

with f s being of the form f s = f s (E S} - 0*sj m s > (VO > A fe i) . The notation here is similar to that adopted in the previous 
section, with h 1 Ds representing the diamagnetic part of /* s , while the leading-order distribution f s is now given by 

/.= 



(2Trk B /Msf 2 (T ]ls ) 1/2 T 



M s (v-V s ) 2 , B 
xexp<^ *\ - — -m' - - \, 30 

which is referred to as the bi-Maxwellian KDF. The form of f s has been prescribed in order to allow the existence of: 

1) A finite non-uniform azimuthal flow velocity V s = e^RQs, where fi s is the rotational frequency (the toroidal 
angular velocity). 

2) A finite temperature anisotropy characterized by kB \ T = k B T ± kB% • 

3) A finite toroidal magnetic field, related just to FLR and diamagnetic effects driven by temperature anisotropy. 

4) A non-uniform number density defined as 



r)s exp 



kBT\\ s 



(31) 



with r]s again denoting the pseudo-density and the function X s being prescribed in terms of the particle effective 
electrostatic energy, the centripetal potential and the azimuthal velocity (see the definition in Paper I). 
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5) Separate treatments of the species and sub-species contributions, for which the previous asymptotic ordering is 
assumed to hold. 

Note that, unlike in the open-surface case, the energy-correction contribution h 2 Ds coming from the er^-expansion 
does not appear in the asymptotic solution (|29[) , as this becomes of higher order than h l Ds . 
In this case the structure functions A s are 

A s = T|| S ,Q S ), (32) 

whose physical meaning has been pointed out in the previous section. As before, to leading order in the asymptotic 
expansion (1291) and in the GK expansion, the structure functions prescribe the fluid fields carried by the bi-Maxwellian 
KDF (pO)) . In particular, for closed nested magnetic surfaces and for strongly-magnetized plasmas, the kinetic con- 
straints give: 

A., = A s (VO + O (e s ) . (33) 

Comparison with Eq. ([26|) shows that the energy dependence no longer appears, so that the fluid fields are only V'-flux 
functions. This is because in the present case, to leading order, the effective potential is itself reduced to a flux- 
function, i.e. $ e Jf = e%[t)- This conclusion is in agreement with the small inverse aspect ratio ordering. 
Hence, the functional forms of the leading-order number density, azimuthal flow velocities and temperatures carried 
by the bi-Maxwellian KDF, are uniquely determined in terms of if). 

Concerning the diamagnetic part h 1 Ds , it can easily be shown that it carries the contributions arising from Taylor- 
expanding the particle toroidal canonical momentum. As a consequence, the diamagnetic KDF is a polynomial 
function of the particle velocity, depending linearly on the thermodynamic forces Aj s , with i = 1,4. From Eg. (1331) and 
following the treatments of Papers I and II, the latter are related to the gradients of the structure functions with respect 
to the magnetic flux ip, and are defined as follows: A% s = 9l g£ s , Ai s = 'g? 1 " , A^ s = 9 l g^ s and = The 
A is , with i = 1,4 carry the contributions due to density, temperature, angular velocity and temperature anisotropy 
gradients respectively. This general form of h 1 Ds follows from the assumed kinetic equilibrium defined by Eg. (12^1) . 
The diamagnetic effects carried by h 1 Ds determine the first-order FLR corrections to the fluid moments of /* s . They 
are important for characterizing collisionless plasmas in the presence of non-uniform differential rotation, density and 
temperature gradients and non-uniform temperature anisotropy. 



VI. THE QUASI-STATIONARY KDF FOR WEAKLY-MAGNETIZED PLASMAS 

In this section we derive a solution for the quasi-stationary KDF describing weakly-magnetized collisionless AD 
plasmas. The basic difference from the strongly-magnetized regime is that gyrokinetic theory does not hold for 
weakly-magnetized plasmas because single particles are effectively not magnetically confined (the Larmor radius is of 
the same order as the characteristic equilibrium scale- length of the plasma, or larger). Hence, contrary to the case for 
strongly-magnetized plasmas, guiding-centre adiabatic invariants can no longer be obtained. Here we retain, however, 
a number of physical features relevant for modelling weakly-magnetized plasmas: 

1) Isotropic temperature: for all of the species it is assumed that the temperature is isotropic. 

2) Nested magnetic flux surfaces: the magnetic field is assumed to allow quasi-stationary solutions with magnetic 
flux lines belonging either to closed or open nested magnetic surfaces. 

3) Azimuthal flow velocity: the plasma is characterized by having a primarily differential azimuthal flow velocity, 
whose leading-order expression is V vs = il s R 2 V(p. 

4) Fluid fields: the collisionless plasma is characterized by non-uniform fluid fields, defined in terms of velocity 
moments of the quasi-stationary KDF. 

5) Kinetic constraints: suitable functional dependences must be imposed so as to ensure that the KDF is an 
asymptotic solution of the collisionless Vlasov equation. For weakly-magnetized plasmas, the kinetic constraints are 
found to differ from those considered before for strongly-magnetized plasmas. They include, in particular, constraints 
on the species angular frequency Q, s . 

6) Analytic form: the KDF is required to be a smooth analytic function. 

Given the requirements l)-6), the solution for the KDF cannot, in general, be a Maxwellian. However, it is possible 
to show that they can all be satisfied by a suitably-generalized Maxwellian, the new solution being expressed only in 
terms of the first integral (|T2"j) and the adiabatic invariant (fT3")) . For clarity of notation, in the following we will label 
with f ws the quasi-stationary KDF for weakly-magnetized plasmas. This has the general form 



fws — fws (E s ,p V s, (E s ) , A fe i) , 



(34) 
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with k > 1 and A being identified with a. In this case, the perturbative expansion is carried out only with respect to 
the variable (E s ). In agreement with the above requirements, a particular solution for f ws is given by: 

IWS ~ /o U 7m n3/2^3/2 eX P ) U—f f ' W 

(2Trk B /M s ) TJ S I k b ± W8 j 

which will be referred to as the Generalized Maxwellian KDF for weakly-magnetized plasmas, with E s and p vs being 
defined respectively by Eqs. (fT3| and f| 12|) . In analogy with the previous treatment, we now introduce the following 
structure functions: 

A ws — {Vwsi T ws: $l ws *) . (36) 

In order that Eq. (|35p defines an adiabatic invariant, the following functional dependences must be imposed on the 
structure functions: 

A ws — A ws (E s ) , (37) 

which represent the kinetic constraints for weakly-magnetized plasmas. Note that, for this configuration, only a 
dependence in terms of the particle total energy is retained. Due to the assumption of having a gravitationally-bound 
plasma, this is the only physically admissible choice for the structure functions. 

Using Eqs. (fT2")) - (fl"3")) . an equivalent representation for f ws is provided by the expression 

, nws / M s (v - V ws ) 2 \ 

^ ~ (2,k B /M s f 2 ^ j " ^ B T WS ) • (38) 



Here V,„ s = Q ws R 2 \7ip, while the function n ws is defined as 



k B T w 



(39) 



with 

X ws = M s 1 ^- + ^n ws - Z s e$t ff . (40) 
2 c 

We conclude that the quasi-stationary KDF is an adiabatic invariant because of the kinetic constraints (j3"T| . 



A. Analytical expansion of /, 



Imposing the kinetic constraints introduces an implicit velocity dependence in the structure functions. This can 
be explicitly dealt with by performing an asymptotic analytic expansion of f ws , similar to that made for strongly- 
magnetized plasmas. In view of the form of Eq. (|37p . a convenient Taylor expansion for f ws can be obtained here 
in terms of the parameter a s ; there is no expansion in the ratio For gravitationally bound plasmas, we can 

set E s — Z s e§ e Jf + O (cr s ) to leading order. Correspondingly, the linear approximation for the structure functions, 
obtained neglecting corrections of O (crj) with k > 2, is 

A ws S A s + (E s - Z 8 e$ e /!) 

where 

A s = [A ws ] Es=Zse<fJ . (42) 
Then, the following relation for f ws holds, correct to first-order in a s : 

Us S /I [1 + h ws ] , (43) 



dA, 



dE s 



(41) 
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where /° s ( y E s ,p ips , (Z s e<& e Jf) , X k t) is the leading-order solution and h ws represents the first-order perturbative con- 
tribution coming from the Taylor expansion. The leading-order solution f^ s can be expressed as 



J WS 



{2vk B /M s ) 6lz T, 



3/2 rpZ/2 



exp 



Ms (v - V 8 
2k B T s 



(44) 



which we will refer to as the drifted Maxwellian KDF for weakly-magnetized plasmas (i.e., a Maxwellian in the species 
co- moving frame having velocity V s ), and which depends on the number density n St temperature T s and azimuthal 
flow velocity V s = ri s i? 2 V(^, with fl s representing the leading-order azimuthal rotational frequency. The leading-order 
number density is then defined as 



n s = rj s exp 



k B T s 



with 



X s = M s 



IV, 



-vn s - Z s e$ 



eff 



(45) 



(46) 



From Eq. (j42|) . to leading-order the structure functions A s = (r) s ,T S) O s ) must satisfy the following kinetic constraints: 

A s = A s (&Jf) + O (a s ) . (47) 

Finally, the first-order correction h ws is again a polynomial function of the particle velocity and has the form: 

Ks = hl s v 2 + h 2 ws v 2 (v • e v ) , (48) 

where we have explicitly singled out the dependences on the particle velocity. The contributions appearing in Eq. (|48[) 
are defined as follows: 



2Z s e 
M 2 Rn, 



ck B T, 



Z s eihVl s j-i 
^3s- 



k B T 3 



D 



2 s 



[D. 



3s 



D 2 



(49) 
(50) 



ws ~ 2Z s ek B T ! 

where we have introduced the following definitions for the energy gradients of the structure functions: D- = 31n a 



'Is 

Z?2s = d Q^} i D^s = d g ^jf ■ These quantities can again be interpreted as generalized thermodynamic forces. 
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B. Properties and discussion 

A number of comments should be made about the properties of the kinetic solution obtained here and its physical 
meaning: 

1) The existence of Eq. ([38j) demonstrates that quasi-stationary drifted Maxwellian kinetic solutions exist also for 
weakly-magnetized collisionless AD plasmas. 

2) Similarly to the case of strongly-magnetized plasmas, the Maxwellian KDF obtained here is generally not an 
exact solution in a strict sense. However, within the validity of the asymptotic expansion (|43j) . this becomes an 
asymptotic equilibrium solution. 

3) The leading-order expressions for the number density, temperature and azimuthal flow velocity appearing in the 
Maxwellian KDF (|44|) . arc found to be functions of 

4) The first-order perturbation h ws allows one to include kinetic effects arising in weakly-magnetized plasmas. 
These contributions are due to the thermodynamic forces D isi i = 1,3, which carry information about the energy 
gradients of the structure functions. 

5) In this case the equilibrium is compatible only with an azimuthal flow velocity, so that accretion flows can only 
occur as a result of turbulence phenomena and cannot be described as equilibrium solutions (see Sections 9 and 10). 

6) The choice (|37[) adopted here for the kinetic constraints is an intrinsic feature of gravitationally bound weakly- 
magnetized plasmas. For this plasma regime, the particle dynamics is mainly determined by the gravitational potential 
and, in principle, also by the electrostatic potential, as shown by Eq. (147p . In contrast, for strongly- magnetized plasmas 
with closed surfaces, the dynamics is determined mainly by the magnetic field and the structure functions are functions 
of the poloidal magnetic flux ip. 
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VII. NUMBER DENSITY AND FLOW VELOCITY 



In this section we address the calculation of the relevant fluid fields associated with the quasi-stationary KDF for 
both strongly and weakly- magnetized plasmas. By definition, given a distribution function /„, a generic fluid field is 
expressed as an integral of the distribution over the velocity space, having the form 



L 



d 3 ^(x)/ s , (51) 



where Z (x) is an arbitrary velocity- weight function and r„ denotes the appropriate velocity space for the integration. 
In this section we focus attention on computing the species number density and flow velocity for plasmas in the high 
and low magnetic field regimes. Using Z (x) = 1 and Z (x) = v, one obtains: 
a) species number density 



n tot - 



= d 3 vf s ; (52) 



b) species flow velocity 



Knowledge of these two fluid moments is required in order to write the Poisson and Ampere equations for studying the 
self-generated EM fields. A basic feature of the present calculation is that the fluid fields arc computed analytically 
in closed form by adopting the asymptotic analytic expansions of the quasi-stationary KDFs for /* s and f ws . Note 
that these velocity moments are uniquely determined once the quasi-stationary KDFs /* s and f ws are prescribed in 
terms of the structure functions A s and A„, s . On the other hand, the equilibrium fluid moments which follow from this 
calculation are identically solutions of the corresponding fluid moment equations. These can be obtained as velocity 
integrals of the Vlasov equation, of the form 



L 



d 3 vZ(x)±f s = 0. (54) 
p clt 



Finally, concerning the notation adopted, here and in the rest of the paper the suffix "tot" is used to label fluid fields 
expressed in terms of f ts or f ws , and to distinguish them from their leading-order solutions computed by means of 
the asymptotic expansions of the same KDFs. 



A. Strongly-magnetized plasmas 

For strongly-magnetized plasmas, the fluid fields must be computed by first performing a transformation of all of 
the guiding-centre quantities appearing in the quasi-stationary KDF to the actual particle position (guiding-centre 
back-transformation). The order of accuracy of this transformation is measured in terms of the parameter em,s 
and depends on the corresponding order of accuracy of the adiabatic invariants used in the solution for the KDF. 
Contributions coming from this transformation carry FLR corrections to the fluid fields, which operate together with 
the FLR-diamagnetic and energy-correction contributions carried by the first-order perturbations of the KDFs. Here a 
distinction must be made between the cases of open and closed magnetic surfaces. In fact, for open-field configurations 
it is possible to show that £m,s J$ while for closed-field configurations, within the validity of the inverse aspect 
ratio ordering, it follows that ■ <~ O (5) <C 1. Hence, FLR effects from the guiding-centre back-transformation are 
negligible for strongly-magnetized plasmas with closed magnetic surfaces. Also, as indicated above, for closed-field 
configurations the terms contributing to the KDF coming from the cr s -expansion are also negligible with respect to 
those proportional to e s and compared with the open-field case. From these considerations it follows that, when 
only first-order contributions are retained in the asymptotic expansion, then the corresponding first-order corrections 
coming from the guiding-centre back-transformation need to be retained only for the leading-order KDF. Finally, 
because of the existence of multiple-species plasmas, which for strong magnetic fields may include also velocity-space 
sub-species, the velocity sub-space T u of integration must be properly prescribed. In fact, charged particles in both 
open and closed configurations can have mirror points (TPs and BPs) or be PPs, which are free to stream through 
the boundaries of the domain. These populations give different contributions to the relevant fluid fields and therefore 
require separate statistical treatments. 
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We consider first the species number density. For the general case of open magnetic surfaces the analytical expansion 
for the quasi-stationary KDF is given by Eq. (|20l) . Then, to first-order in all of the expansion parameters (e s , a s and 
Em,s), the number density is given by 

nt ot = / d 3 v{[f s ] GK [1 +h 1 Ds + h 2 Ds ]}, (55) 
where f s denotes the leading-order KDF to which the guiding-centre back-transformation must be applied up to 

L J GK 

first-order in Sm.s- The corresponding expression holding for closed surfaces in the case of small inverse aspect ratio 
ordering reduces to 



nt ot 



jf d 3 v{f s [l + h 1 Ds ]}, (56) 



where we have made use of Eq. (|29|) . In both cases the constitutive equation for the total number density can be 
written as 

n*°* a n s [1 + A„J , (57) 

where the leading-order contribution n s is given by Eq. (|3"Tj) for closed surfaces and by Eq. (f22"]) for open surfaces when 
the ordering Em.s *C e s holds, while A„ s carries the first-order corrections associated with FLR, diamagnetic and 
energy-correction effects. 

Next, let us consider the species flow velocity. Velocity-space integrals analogous to (|55|) and ([56]) can be written 
also in this case. The species flow velocity can be generally represented in terms of the constitutive equation 

nl ot Vl ot = n s [U s + A Us ] . (58) 

For open magnetic surfaces and when em,s *C e s , U s = V s + U\\ s h is the leading-order flow velocity carried by the 
bi-Maxwellian KDF with a parallel velocity perturbation, in which V s = tt s R 2 Vtp and f7ii s = -g£ s - Moreover, the 
two frequencies Cl s and £ s are subject to the kinetic constraints given by Eq. (|26[) . The second term Au s represents 
the first-order correction which can be decomposed as follows: 

A Us = (A UsV , W x Vip, A Us¥ ,V^, A Usl5 V?/> x V<^) . (59) 

Note that for open magnetic surfaces: 

1) All components of Au, are linear functions of the thermodynamic forces appearing in the first-order perturbation 
of the KDF. 

2) The component Au s¥ > provides a correction to the leading-order azimuthal velocity. This is non- vanishing even 
in the absence of any guiding-centre contribution in the KDF and also in the case of isotropic temperature. 

3) The component Au s i? is related to the temperature anisotropy as well as to both FLR-diamagnetic and energy- 
correction effects. 

4) The component Au s i /, is associated with FLR effects coming from the guiding-centre back-transformation. 

5) Both J7ii s b and Au s can give rise to inward or outward flows of matter, both in the radial and vertical directions. 
The existence of a non- vanishing parallel velocity U\\ s in the kinetic solution is allowed by the gyrokinetic conservation 
laws. This means that it can always be suitably prescribed in agreement with the kinetic constraints. On the other 
hand, the correction Au, acquires a precise physical meaning within the present formulation. In fact Au s is generated 
by the existence of a non-uniform and non-isotropic plasma. Its precise form is determined automatically once the 
kinetic constraints for the leading-order structure functions are prescribed. 

6) The expression for Au, depends on the particle sub-species (i.e., the distinction between PPs, TPs, and BPs). In 
fact, each of these populations gives different contributions to the components of the velocity Au, . As indicated below, 
they give rise to interesting physical phenomena for the disc dynamics, in relationship with the Ampere equation and 
the self-generation of both toroidal and poloidal equilibrium magnetic fields. 

We consider now the case of closed magnetic surfaces. The calculation is here made simpler because, by assumption, 
U\\ a vanishes and only the e s -expansion is relevant in the asymptotic expansion. Thus, the flow velocity can still be 
written as in Eq. (|58[) . but now only an azimuthal leading-order flow velocity U s = f2 s i? 2 V</? can arise, with fl s = fl s 
as indicated above. On the other hand, the first-order correction Au s now reduces to 

A Us = (0, A Us¥ ,V(^, Au^Vip x Vip) . (60) 
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In particular, the V$ x component vanishes in this approximation since it is related to terms coming from the 
guiding-centre back transformation which are of higher-order for closed surfaces. Again, the poloidal component Au st j 
is due to diamagnetic FLR velocity corrections produced by temperature anisotropy (see also Paper I). However, under 
the hypothesis of closed nested magnetic surfaces, it cannot give rise to a net accretion velocity. 

We therefore conclude that in both cases plasma temperature anisotropy affects the existence of non-vanishing 
species poloidal flow velocities. This feature only occurs for strongly-magnetized plasmas. The importance of this 
result for AD plasmas lies in the fact that equilibrium poloidal flow velocities may also give rise to a net poloidal current 
density. The latter in turn will generate a finite equilibrium toroidal magnetic field. Therefore, species temperature 
anisotropies in collisionless AD plasmas actually provide an effective physical mechanism for the self-generation of 
toroidal magnetic field in these systems. 



B. Weakly- magnetized plasmas 

For weakly-magnetized plasmas, the explicit calculation of the fluid fields is made simpler by the fact that the 
kinetic equilibrium does not contain any guiding-centre adiabatic invariant. In this case, only contributions arising 
from the o^-expansion need to be taken into account. To first-order in a s the species number density can then be 
easily calculated as follows: 



Sv{f ws [l + h ws ]}, 



(61) 



which recovers again the constitutive equation 



= 71, fl + A, 



(62) 



Here the leading-order term n s is defined by Eq. (|45| and is the contribution carried by the Maxwellian KDF. On the 
other hand, the first-order correction A„ s is found to be given by 



A, 



V? 



3k B T s 



9k B T s 
M, 



(63) 



where V 2 = Q 2 R 2 and the quantities h^ s and h 2 us are defined by Eqs. (|49)) - ((50j) . Similarly, in the same approximation, 
the flow velocity is given by the velocity integral 



n s JT„. 



(64) 



After explicit calculation, this gives 



n tot Y tot 



n s [V s 



(65) 



where V s = il s R 2 Vip is the leading-order flow velocity carried by the Maxwellian KDF. The first-order contribution 
Av s is similarly found to be given by 



V s A„ s 



+K 



!5k B T s 
M, 



Qk B T s 
M s 

f V 2 + 2Q 2 R'' 



k B T s 



(66) 



The following features should be noted: 

1) The flow velocity is purely azimuthal. No further components of the flow velocity are allowed for weakly- 
magnetized plasmas, in contrast with the case of strongly-magnetized plasmas. 

2) When the gravitational potential energy is dominant over the electrostatic energy, the functional dependence of 
the structure functions as implied by the kinetic constraints is determined primarily by the gravitational potential. 
For example, the leading-order azimuthal rotational frequency Q s is simply of the form n s ~ Q s (^g)- This means 
that its functional dependence is compatible with the Keplerian rotational frequency. 
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VIII. PRESSURE TENSOR AND EQUATIONS OF STATE 



In this section we focus on the calculation of the pressure tensor corresponding to the kinetic equilibria obtained 
for both strongly and weakly-magnetized plasmas. The species pressure tensor (or partial pressure tensor) is defined 
with respect to the species flow velocity as the velocity moment 

n tot = / d 3 vM s (v - V* ot ) (v - V* ot ) (67) 



In the context of the present kinetic treatment n tot is then uniquely prescribed in terms of the quasi-stationary KDF. 
As indicated below, this enables us to determine to the requisite accuracy also the corresponding equations of state 
relating the components of the tensor to the structure functions, giving a self-consistent treatment of the physical 
properties of quasi-stationary collisionless AD plasmas. In the following we provide explicit expressions for n tot ; the 
overall pressure tensor of the system can then be obtained by summing over the contributions from the separate 
species: g* * = E s=i , e flf • 

We consider first the case of strongly-magnetized plasmas, adopting respectively for configurations with open and 
closed magnetic surfaces the asymptotic expansions given by Eqs. ([2T)l) and (j2"9"j) . As a result of the perturbative 
calculation, the species pressure tensor is represented as 

Sf* -S s + A E S > ( 68 ) 

where II is the leading-order term (with respect to all of the expansion parameters), while An represents the 
first-order correction. For both closed and open magnetic surfaces the tensor II is obtained using a bi-Maxwellian 
KDF with temperature anisotropy. To represent n , we introduce for convenience the set of right-handed orthogonal 

unit vectors (b,ei,e2), where b = ||, while ei and e 2 are two orthogonal vectors in the plane perpendicular to 
the magnetic field. In terms of this basis, the unit tensor can be represented as: I = bb + + e 2 e 2 . Then, it 
follows that with respect to the unit tensor I, the pressure tensor II is symmetric, "diagonal and non-isotropic, with 
a representation of the form: 

|I s = n s k B T Xs l + n s k B (T\ [s - T Xs ) bb. (69) 

Here n s is the leading-order species number density and Xju, T± s are the leading-order species parallel and perpen- 
dicular temperatures. Instead, the precise form of An is geometry-dependent and contains FLR corrections. It is 
always possible to represent it in terms of the general decomposition: 

Ag s = AI^I + AIT>b + A1J 3 , (70) 

in which An* and AII^ are diagonal first-order anisotropic corrections to the pressure tensor, while An 3 in this basis 
is generally non-diagonal. For st rong ly- magnetized plasmas and closed magnetic surfaces, the precise form of the 
tensor pressure has been given in [22j j . The physical properties of the solutions (|69[) and (|70[) can be summarized as 
follows: 

1) The total tensor pressure II is symmetric. 

2) The leading-order pressure tensor II calculated in this approximation is diagonal but non-isotropic. We note 
that the source of anisotropy in Eq. (|6"9")) is provided by the temperature anisotropy. 

3) The first-order correction An instead is generally non-diagonal and non-isotropic in the (b, ei,e2) basis. Two 
different physical mechanisms contribute to generating this effect. The first one is again the temperature anisotropy, 
while the second is produced by first-order perturbative corrections to the KDF, and so depends linearly on the 
thermodynamic forces. 

We next consider the case of weakly-magnetized plasmas. Performing a similar calculation, it is possible to prove 
that, to first-order in cr s , the pressure tensor is symmetric and isotropic and can be written as 

2* ot ~ p s +An s ]|. (71) 

In particular, the leading-order term is defined as 

II S = n s k B T s , (72) 



with the number density being given by Eq. pS)) . Both the number density n s and the temperature T s are subject to 
the kinetic constraints expressed by Eg. (1471) . The explicit representation of the first-order isotropic term AII S is as 
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follows: 

15fc|T s 2 2 3k B r s - 
M1 s M s 

+h 2 ws n s ^lsR 

A striking feature of Eq. (|73[) is the explicit dependence in terms of the thermodynamic forces, which take into account 
the gradients of the structure functions. Again, it should be noted that the pressure tensor (|7ip is isotropic because 
the associated KDF has an isotropic temperature and does not contain guiding-centre or FLR effects. 



An s = h\,„n a 



M? 



V. 



2 Sk B T s 
Ms 



(73) 



A. Equilibrium equations of state 

An important result of the present theory is the explicit construction of equations of state for the various components 
of the species pressure tensor. For definiteness, let us consider here only the leading-order contributions to with 
respect to the relevant expansion parameters. This provides finite-term equations (leading-order equations of state) of 
the form n*°* = n' ot (A s ), with A s being the appropriate structure functions. Note that in principle the solution for 

n*°* allows one to obtain equations of state which include also first-order corrections which are linearly proportional 
to the thermodynamic forces. 

In particular, in the case of strongly-magnetized plasmas this recovers for the leading-order perpendicular and 
parallel pressures the expressions 

p ±s = n s k B T ±s , (74) 
p\\ s ee n s k B T\\s- (75) 

The number density and the parallel and perpendicular temperatures are subject to their respective kinetic constraints 
(see the discussion above). Eqs. (j22|) and (|3Tj) provide a clear representation of the physical effects contributing to 
the equations of state for strongly-magnetized plasmas in the case of open and closed magnetic surfaces respectively. 
In particular, the functions X s allow corrections due to particle electrostatic and gravitational energy, centripetal 
potential and azimuthal and parallel flows to be explicitly taken into account. 

An analogous equation of state can be obtained in the case of weakly-magnetized plasmas for the isotropic pressure 
tensor: 

p s ee n s k B T s , (76) 

with p s denoting the species scalar pressure and n s being given by Eq. (|45p to leading-order. This equation of state 
allows one to clearly display the contributions due to the gravitational and electrostatic potentials as well as the 
azimuthal flow velocity. 

Note that Eqs.([74")) and ([75]). as well as Eq. (jT6"l) . provide only the leading-order solution for the corresponding 
equations of state. In fact, a more accurate solution should necessarily include also higher-order terms coming from 
the diamagnetic and energy-correction contributions, both of which result from the Taylor expansions performed on 
the KDF. These corrections are responsible for the appearance of the terms An and AIT in Eqs. ([68|) and ([7T]) 
respectively. This implies that, for both strongly and weakly magnetized collisionless plasmas, the equations of state 
for the partial pressures cannot be of the frequently-used polytropic type p — Kp 1 with k and 7 both being constants 
(and p again here being the mass density). This is not surprising because the quasi-stationary kinetic solutions 
considered here are not thermodynamic equilibria (as required for deriving this polytropic equation of state from 
microscopic considerations). This is clear from the appearance of non- vanishing thermodynamic forces. 

The present kinetic approach leads naturally to the use of temperature and density as thermodynamic variables, 
since these fluid fields are directly related to the (leading-order) structure functions contained in the equilibrium KDF. 
On the other hand, for a proton-electron plasma, which does not possess internal binding energy, the temperature 
also represents a statistical measure of the specific internal energy of the system, which by definition does not include 
potential or EM energy. We note that equivalent representations can be given in terms of temperature or of specific 
internal energy although the more general case with a non-isotropic pressure tensor requires the introduction of the 
concept of "directional" specific internal energy (related to the anisotropic temperatures). 
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IX. ANGULAR MOMENTUM 



In this section we discuss the implications of the kinetic treatment for the law of conservation of fluid angular 
momentum. For doing this, we first define the species fluid canonical toroidal momentum as 



1 

■*tot 



d 3 



Z s e 
c 



(77) 



for f s — /* s or f s — f ws respectively in the cases of strongly and weakly-magnetized plasmas. Consider then the 
corresponding conservation law for the species total canonical momentum. This can be recovered by identifying the 
weight function Z — ip rs in Eq. (j54")l . i.e. by setting 



j dt 



0. 



In the equilibrium case this implies the species fluid angular momentum conservation law 

Z^e. 



^tot\rtot t tot 
l a V D 



-V^ • nf'Vf * = 



for the species angular momentum 



(78) 



(79) 



(80) 



where expressions for the number density, flow velocity and pressure tensor have been derived in the previous sections. 
In Eq. (f79"]) a key role is played by the divergence of the species pressure tensor. For strongly- magnetized plasmas, 
using the leading-order expression (f69|). this is given by: 



bB V 



P\\, 



PA 



B 



- A Ps Q, 



(81) 



where Q ee [bb • Vlni? + |g-b x J — VlnS] and Ap s ee (pu s — p± s )- It is clear that in this case V • n tot has non- 
vanishing components in arbitrary spatial directions, including the azimuthal direction along V</?. On the other hand, 
for weakly-magnetized plasmas, Eq. ([7T|) gives 



v • n tot = v [n s + An s ] ■ i. 



(82) 



Since the pressure tensor is isotropic and we are assuming axisymmetry, it follows that the component of V ■ II ° 
along V</? must vanish identically. 

For a single species, the total canonical momentum L*°* and the total angular momentum L* ot in general differ 
because of the contribution of the magnetic part proportional to the flux function tp. However, a different conclusion 
can be drawn for the corresponding canonical momentum density n*°*L*°* and angular momentum density n*°*L' ot . 
If one considers summation over species for both these quantities and imposes the quasi-neutrality condition 



0, 



(83) 



then one obtains the identity 



(84) 



We next investigate the consequences of Eq. (l79)) for the dynamical properties of collisionless plasmas. Note the 
following aspects: 

1) In the usual interpretation and also for weakly- magnetized collisionless plasmas within our treatment (see Eq. (|85[) 
below), the directional derivative of L* ot along the flow velocity V' ot vanishes. However, for strongly-magnetized 
plasmas Eq. (|79[) shows that equilibrium configurations are possible in which this is generally non-zero. This arises 
because of the non-isotropic pressure tensor and the poloidal components of the flow velocity which, in turn, are 
consequences of temperature anisotropy, the first-order energy-correction and FLR-diamagnetic effects which are not 
included in standard MHD treatments. 
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2) According to Eq. (|79p . spatial variation in the species angular momentum implies the possibility of having quasi- 
stationary radial matter flows in the disc without departing from the unperturbed equilibrium solution. These can 
correspond either to local outflows or inflows; both can be described consistently within the present kinetic solution for 
open magnetic field lines in strongly-magnetized plasmas. Local inflows and outflows can occur independently and are 
described consistently by their respective quasi-stationary KDFs. Radial flows arise due both to the parallel velocities 
Uu s and to the kinetic effects driven by the first-order energy-correction and FLR-diamagnetic effects. Therefore, 
species radial flows appear necessarily together with a non-isotropic pressure tensor and a non-vanishing toroidal 
magnetic field (see also Papers I and II and the discussion below). 

3) For weakly-magnetized plasmas, the tensor pressure is isotropic and Eq. (|79[) reduces to 

n totytot . VL tot = Q) (g 5 ) 

which, thanks to axisymmetry and recalling Eq. (|65|) . is identically satisfied. Hence, in this case equilibrium radial 
flows are excluded. 

In conclusion, for the case of weakly- magnetized plasmas, it follows from the present treatment that quasi-stationary 
equilibrium configurations with isotropic pressure tensor would not have any net radial flow. This is as expected. In 
order to have an accretion flow in this context one needs to have some form of effective viscosity, either appearing 
explicitly or coming from perturbations around the equilibrium state such as those leading to MRI in the conventional 
picture. A particular goal of our future work will be to investigate perturbations around the equilibrium states 
presented here, to see whether a process analogous to MRI then appears also for collisionless plasmas. 

On the other hand, in the case of strongly-magnetized plasmas the situation is different and net radial inflow can 
occur even in the absence of effective viscosity being explicitly added or coming from perturbative mechanisms. For 
strongly-magnetized plasmas, Eq. (|79[) implies that particles of one species can move radially in a quasi-stationary 
configuration independently of those of other species, with the species angular momentum not being conserved, since 
the conservation law involves the canonical momentum (including a magnetic-field contribution) and not just the 
standard angular momentum. The flow velocity V*°* is different for each species and, also, the pressure tensor is non- 
isotropic. Hence, even with quasi-neutrality, the total angular momentum of the matter can change due to balance 
with the torque produced by the non-isotropic species pressure tensor and/or a net current flow across magnetic 
surfaces (see also the discussion below) . As indicated above, the non-isotropic nature of n tot is caused by temperature 
anisotropy as well as by first-order and FLR effects, which are determined self-consistentfy in the kinetic approach. 



X. THE KINETIC ACCRETION LAW 



In this section we discuss how quasi-stationary accretion flows could occur in collisionless AD plasmas as a result 
just of the equilibrium configuration without requiring additional effective viscosity of the sort mentioned above. Here 
the role of viscous stresses is played by the anisotropic pressure tensor, which is part of the equilibrium solution. 

For the physical conditions considered in the Introduction to which the theory applies, the characteristic time for 
the inward accretion flow in accretion discs is typically longer than the characteristic Larmor time tl s as well as the 
Langmuir time r ps and smaller than the Spitzer ion collision time so that an accretion flow can be consistently described 
by the present collisionless kinetic treatment for quasi-stationary equilibria. In this section we will demonstrate that 
equilibrium accretion flows cannot arise in the case of strongly-magnetized plasmas with closed magnetic surfaces and 
weakly-magnetized plasmas. Instead, equilibrium accretion flows are permitted for general open-field configurations. 
Note, however, that in general there are both open and closed surfaces co-existing, as illustrated in Fig.l. 

Let us consider first domains which are locally characterized by open flux surfaces. In such domains, parallel flows 
can be included in the quasi-stationary KDF only if the guiding-centre canonical momentum p' vs is conserved, i.e. the 
plasma is strongly-magnetized (see Section 4). In fact, in such cases, the quasi-stationary KDF (|20p can sustain both 
poloidal and radial species flow velocities, which are defined respectively as 

V ps = V* ot • = -L [ d 3 v[v e„] £, (86) 

V Rs ee vr-e fl = -^4t (87) 
where e# = y^gy, e R = and the species mass radial current density J R °* is defined as 

Jrs=J t d 3 v[ve R ]fc s . (88) 
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In the velocity-space integrals indicated above, the contributions from PPs, BPs and TPs need to be distinguished. 

The physically-relevant situations are those in which there is a non-vanishing net radial species accretion flow, i.e. 
where the average species radial mass current ((Jr s )) = z l z J z ^ Jnsdz is negative (inward flow), with z\ and Z2 
being suitably prescribed. Although local contributions to ((Jr s )) can arise from TPs, BPs and PPs, the overall 
accretion flow is mainly associated with PPs. Notice also that, to leading-order, the presence of poloidal accretion is 
necessarily associated through the expression of Ut\ s to the existence of a toroidal magnetic field. Note the following 
basic features involved in the accretion process. The ratio 4r- is approximately constant due to the kinetic constraints, 
and so 

K.]i ~ [ T p]i 



(89) 



for any two arbitrary positions "1" and "2" prescribed in terms of the magnetic coordinates (tp, Then consider the 
case £m,s •C s s , which allows one to approximate the guiding-centre quantities with the expression for them evaluated 
at the particle position. Assuming that I — I (ip) (see Paper II) and considering the two positions (ip, $i) and (tp, #2) 
on the same flux surface, the kinetic accretion law follows 

In this case, under the same assumptions, it follows from the continuity equation that the ratio of the corresponding 
species number densities must vary on a given ^-surface according to the following relation: 

n s (*M 2 ) B 2 (tp, $2) T\\ s (tp, i?i ) ' 1 ' 

Therefore, on a given ^-surface: 

1) the species parallel flow velocity increases with the parallel temperature while decreasing with respect to the 
magnitude of the magnetic field; 

2) the species number density instead increases with the magnetic pressure and decreases with the parallel temper- 
ature. 

The physical interpretation for both of these is clear: higher magnetic pressure slows down the matter accretion rate 
while increasing the number density, whereas higher parallel temperature corresponds to higher radial fluid mobility, 
thus decreasing the local species number density. 

Next, we consider domains of strongly-magnetized plasmas with closed field lines, again using the inverse aspect 
ratio expansion. We want to prove that the poloidal-averaged radial flow velocity (Vr s )^ vanishes identically, when 
Vr s is computed in terms of the quasi-stationary KDF (f2U|). Given a generic function of the form C (tpj'd), the operator 
(C (ip,!}))^ is defined as 

{c ^^U^W^\ ci ^ h (92) 

with k = J 2?r B rf y^| and so, from Eq. flSfJl) it follows that, to leading-order in the inverse aspect ratio 

which vanishes because, to leading-order in 8, B p - VR = B p Vd-VR is antisymmetric with respect to the transformation 
•d — > -d + 7r. Higher-order corrections in the inverse aspect ratio can be included, but they require also performing the 
asymptotic expansion of the KDF to higher-order in both a s and e s . Therefore, in sub-domains of the plasma where 
magnetic surfaces are closed and nested no net equilibrium accretion flow can arise. 

Finally, the corresponding treatment for weakly-magnetized plasmas can be recovered from Eq. (|43p . Since the 
species flow velocity V* ot has only an azimuthal component, it follows that Vr s is identically zero. Hence, under the 
present assumptions, no net quasi-stationary radial flow can arise in the case of weakly-magnetized plasmas. 

To summarize: the present theory provides a possible new collisionless physical mechanism giving an equilibrium 
accretion process in AD plasmas. In particular, we note that: 

1) Only strongly-magnetized plasmas with open magnetic surfaces can sustain these equilibrium accretion flows. 

2) The primary source of this equilibrium accretion flow mechanism is the appearance of equilibrium radial flows 
driven by temperature anisotropies and phase-space anisotropies. These are directly connected with the existence of 
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non-isotropic species pressure tensors, which in turn play the role of an effective viscosity in driving quasi-stationary 
accretion flows. 

3) Quasi-stationary accretion flows are consistent with the basic conservation laws (for mass density and canonical 
momentum) and with the existence of a non-isotropic species pressure tensor (see also the discussion in Section 11). 

4) The accretion law could, in principle, be tested experimentally if one had suitable observations, since it relates 
the magnitude of the species parallel flow velocity Un a to the local values of the magnetic field magnitude and the 
parallel temperature. 

5) First-order (as well as higher-order) perturbative corrections, can in principle be included consistently in the 
present theory. 



XI. THE KINETIC DYNAMO 



Here we address the problem of the self-generation of magnetic field in quasi-stationary collisionless AD plasmas. 
We refer here to this phenomenon of self-generation of both poloidal and toroidal magnetic fields as a quasi- stationary 
kinetic dynamo effect. This is described by the Ampere equation with current density J defined as 



for f s = f* s or f s — f ws respectively in the cases of strongly and weakly-magnetized plasmas. In all cases considered 
here, the current density can be generally represented in terms of the magnetic coordinates {tp 1 (p 1 , d) as 

J = ( J V ,V?9 x Vip, J v Vip, J#Vip x Vip) . (95) 

In particular, based on the calculation of the flow velocity (see Section 10), it can be shown that: 

a) For strongly-magnetized plasmas with open magnetic surfaces all of the three components (J^,, J v , </$) are 
generally non- vanishing. These include the contributions carried by TPs, BPs and PPs. 

b) For strongly-magnetized plasmas with closed magnetic surfaces, vanishes identically. In this case the current 
includes contributions from both TPs and PPs. 

c) For weakly-magnetized plasmas both J.^ and J$ vanish identically. 

The toroidal component of the Ampere equation gives the generalized Grad-Shafranov equation for the poloidal 
flux function 

4-7T 

A*V P = - — J v , (96) 

where A* = R 2 V ■ (i?~ 2 V) is an elliptic differential operator. The components of Ampere's equation along the 
directions Vi9 x Vip and Vip x \7ip provide instead two PDEs for the toroidal component of the magnetic field 
|B T | = I/R (see the definition in Eq.0): 

£ = W 
oip c 

dl 4tt 

M = (98) 



which give the general solubility constraint 



dJtf) dJ{) 



(99) 



which is equivalent to imposing the charge continuity equation V • J = 0. In the case of closed surfaces, the additional 
solubility condition 

J^dd = (100) 



must also be imposed. Both constraint equations must be taken as being solubility conditions to be satisfied due to 
the arbitrariness of the structure functions. 
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For open surfaces, the function / is of the general form I(ip,d) and so Eq. (|99|) can always be satisfied. Instead, 
for closed surfaces / = / (ijj), which can always be asymptotically satisfied if inverse aspect ratio ordering applies. 
Constraints ([99]) and (jlOOp are then both satisfied. 

Various scenarios can be envisaged in which quasi-stationary kinetic dynamos can be present. We will list the basic 
features of this mechanism: 

1) In all of the cases discussed above, J v is generally non- vanishing, implying the existence of a self-generated 
poloidal magnetic field. Instead, a toroidal component of the magnetic field only arises if a poloidal current is present. 

2) For strongly-magnetized plasmas, first-order FLR-diamagnetic and energy-correction effects, driven by tem- 
perature anisotropy, are responsible for the generation of poloidal currents, and hence of toroidal magnetic field. 
Gyrophase-dependent contributions can arise in this case, driven by the same thermodynamic forces. These originate 
from the guiding-centre back-transformation, which is characteristic of open-field configurations, and are responsible 
for the generation of J^. For strongly-magnetized plasmas with open magnetic surfaces, the parallel velocity U\\ s also 
contributes to the generation of poloidal currents. 

3) For weakly-magnetized plasmas, only the poloidal component of the magnetic field can be self-generated. How- 
ever, when the azimuthal angular velocity coincides to leading-order with the Keplerian frequency, the azimuthal 
current density vanishes due to quasi- neutrality (to leading-order). Therefore, in this case the current is necessarily 
produced by first-order corrections. This conclusion is consistent with the assumption of weak magnetic field. 

We stress that, in contrast to customary MHD treatments, the quasi-stationary kinetic dynamo effect described 
here can occur even in the absence of possible instabilities or turbulence phenomena. In particular, configurations 
with closed magnetic surfaces or contributions from TPs in the case of open surfaces could be responsible for the 
self-generation of toroidal field even without needing any net accretion flow in the domain of interest. This toroidal 
field is associated with the existence of kinetic torques which cause redistribution of angular momentum, as discussed 
above. 



XII. CONCLUSIONS 

In this paper a kinetic description for collisionless quasi-stationary accretion disc plasmas has been formulated within 
the framework of Vlasov-Maxwell theory. A perturbative approach has been developed, which enables the systematic 
analytical construction of equilibrium kinetic distribution functions for non-relativistic axisymmetric collisionless 
plasmas subject to both gravitational and electromagnetic fields. The cases of both weakly and strongly-magnetized 
plasmas have been investigated, for configurations with both open and closed magnetic surfaces. The main aim was to 
establish kinetic equilibrium configurations to use as a starting point for subsequent perturbation analysis, looking for 
kinetic mechanisms which could give rise to an effective viscosity able to drive accretion flows within the collisionless 
regime. 

Equilibrium solutions for the Vlasov equation with non-uniform number density, azimuthal rotation, possible accre- 
tion flows and non-uniform temperature anisotropy have been constructed. For doing this, the distribution function 
has been represented in terms of first integrals and adiabatic invariants, as follows from conservation laws of the single- 
particle dynamics, and taking into account suitable kinetic constraints. As a consequence, based on the perturbative 
kinetic approach, explicit constitutive equations for the fluid fields have been determined, which are accurate to first- 
order in the relevant expansion parameters. This permits the construction of asymptotic solutions of the corresponding 
MHD fluid equations systematically retaining all first-order corrections, including effects due to FLR-diamagnetic and 
energy-correction contributions. 

Several physical issues have been analyzed. These concern: a) deriving equilibrium equations of state for the species 
pressure tensor components; b) establishing a fluid angular momentum conservation law and comparing it with the 
predictions of standard fluid treatments; c) investigating possible kinetic accretion processes within equilibrium AD 
configurations (without perturbing the equilibrium) and deriving a related kinetic accretion law; d) demonstrating the 
existence of a quasi-stationary kinetic dynamo mechanism for the self-generation of poloidal and toroidal magnetic 
field. 

Potential applications of the theory include, in particular, the case of radiatively inefficient accretion flows. As 
mentioned above, future work will address making stability analysis of the kinetic equilibria presented here, so as to 
investigate further mechanisms which could be responsible for driving the accretion flow in the collisionless regime. 

Acknowledgments - This work has been partly developed in the framework of MIUR (Italian Ministry for 
Universities and Research) PRIN Research Programmes and of the Consortium for Magnetofluid Dynamics, Trieste, 
Italy. 



25 



[1] Narayan R., Mahadevan R., Quataert E., 1998, in Theory of Black Hole Accretion Discs, 148, ed. M. Abramowicz, G. 

Bjornsson and J. Pringle, Cambridge University Press, Cambridge (UK). 
[2] Saxton C.J., Wu K., Cropper M., Ramsay G., 2005, Mon. Not. R. Astron. Soc., 360, 1091-1104. 
[3] Saxton C.J., Wu K., Canalle J.B.G., Cropper M., Ramsay G., 2007, Mon. Not. R. Astron. Soc., 379, 779-790. 
[4] Cremaschini C, Miller J.C., Tessarotto M., 2010, Phys. Plasmas, 17, 072902. 
[5] Cremaschini C, Miller J.C., Tessarotto M., 2011, Phys. Plasmas, 18, 062901. 

[6] Frank J., King A., Raine D., 2002, Accretion power in astrophysics, Cambridge University Press, Cambridge (UK). 

[7] Vietri M., 2008, Foundations of High-Energy Astrophysics, University Of Chicago Press, Chicago, USA. 

[8] Coppi B., 2005, Phys. Plasmas, 12, 057302. 

[9] Coppi B., Rousseau F., 2006, Astrophys. J., 641, 458-470. 
[10] Cremaschini C, Beklemishev A., Miller J., Tessarotto M., 2008, AIP Conf. Proa, 1084, 1067-1072. 
[11] Cremaschini C, Beklemishev A., Miller J., Tessarotto M., 2008, AIP Conf. Proa, 1084, 1073-1078. 
[12] Kocharovsky V.V., Kocharovsky Vl.V., Martyanov V.Ju., 2010, Phys. Rev. Letters, 104, 215002. 
[13] Catto P.J., Bernstein I.B., Tessarotto M., 1987, Phys. Fluids B, 30, 2784. 

[14] Cremaschini C, Miller J.C., Tessarotto M., 2010, Theory of quasi- stationary kinetic dynamos in magnetized accretion 
disks, Proceedings of the International Astronomical Union No. 274, Advances in Plasma Astrophysics, Giardini Naxos, 
Sicily, Italy, Sept. 06-10, 2010, Cambridge University Press (Cambridge, U.K.), vol. 6, p. 228-231. 

[15] Blandford R.D., Payne D.G., 1982, Mon. Not. R. Astron. Soc., 199, 883. 

[16] Rebusco P., Umurhan O.M., Kluzniak W., Regev O., 2009, Phys. Fluids, 21, 076601. 

[17] Mikhailovskii A.B., Lominadze J.G., Churikov A. P., Pustovitov V.D., 2009, Plasma Physics Reports, 35, 4, 273-314. 
[18] Bernstein I.B., Catto P.J., 1985, Phys. Fluids, 28, 1342. 
[19] Littlejohn R.G., 1979, J. Math. Phys., 20, 2445. 
[20] Littlejohn R.G., 1981, Phys. Fluids, 24, 1730. 
[21] Kruskal M., 1962, J. Math. Phys. Sci., 3, 806. 

[22] Cremaschini C, Miller J.C., Tessarotto M., 2010, Kinetic closure conditions for quasi-stationary collisionless axisymmetric 
magnetoplasmas, Proceedings of the International Astronomical Union No. 274, Advances in Plasma Astrophysics, Giardini 
Naxos, Sicily, Italy, Sept. 06-10, 2010, Cambridge University Press (Cambridge, U.K.), vol. 6, p. 236-238. 



